Investigation of heavy-heavy pseudoscalar mesons in thermal QCD Sum Rules 
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We investigate the mass and decay constant of the heavy-heavy pseudoscalar, B c , rj c and r/b 
mesons in the framework of finite temperature QCD sum rules. The annihilation and scattering 
parts of spectral density are calculated in the lowest order of perturbation theory. Taking into 
account the additional operators arising at finite temperature, the nonperturbative corrections are 
also evaluated. The masses and decay constants remain unchanged under T = 100 MeV, but after 
this point, they start to diminish with increasing the temperature. At critical or deconfinement 
temperature, the decay constants reach approximately to 38% of their values in the vacuum, while 
the masses are decreased about 5%, 10% and 2% for B c , r/ c and r\b states, respectively. The results 
at zero temperature are in a good consistency with the existing experimental values as well as 
predictions of the other nonperturbative approaches. 

PACS numbers: 11.55.Hx, 14.40.Pq, ll.10.Wx 
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I. INTRODUCTION 

Over the last two decades, there is an increasing interest on properties of hadrons under extreme conditions [l], Qj. 
According to these investigations, two theoretical aspects, namely theoretical studies of hadrons at finite temperature 
and density as well as a careful analysis of the heavy ion collision results are important. Calculation of hadronic 
parameters at finite temperature and density directly from QCD is a difficult problem. The thermal QCD is successful 
theory in the large momentum transfer region, where the quark-gluon running coupling constant is small and one 
can reliably use perturbative approaches. However, at the hadronic scale, this coupling constant becomes large and 
perturbation theories fail. Hence, investigation of hadronic properties requires some nonperturbative approaches. 
Some nonperturbative approaches are lattice QCD, heavy quark effective theory (HQET), different quark models and 
QCD sum rules. Among these approaches, the QCD sum rule method [3j and its extension to the finite temperature 
[j] has been extensively used as an efficient tool to hadron physics The same as QCD sum rules in vacuum, the 
main idea in thermal QCD sum rules also is to relate the hadronic parameters with the QCD degrees of freedom. In 
this method, an appreciate thermal correlator is expressed in terms of interpolating currents of participating particles. 
From one side, this correlation function is evaluated saturating it by a tower of hadrons with the same quantum 
numbers as the interpolating currents. On the other hand, it is calculated via the operator product expansion 
(OPE) in terms of operators having different mass dimensions. Matching these two different representations of 
the same correlation function provides us a possibility to predict hadronic properties in terms of finite-temperature 
perturbation theory and long-distance nonperturbative physics including the thermal quark and gluon condensates 
as well as thermal average of energy density. 

Comparing to the QCD sum rules in vacuum, the thermal QCD sum rules have several new features. One of them is 
to take into account the interaction of the currents with the existing particles in the medium. Such interactions require 
modification of the hadronic spectral function. The other aspect is breakdown of Lorentz invariance by the choice 
of reference frame. Due to residual 0(3) symmetry at finite temperature, more operators with the same dimensions 
appear in the OPE compared to those at zero temperatu re I6| — 18|] . The thermal QCD sum rule method has been 
extensively used to study the thermal properties of light 0-[Hl|7 heavy-light fl2l - fl4| and heavy- heavy flB - fl8t mesons 
as a reliable and well-established method. 

The discussion of heavy mesons properties at zero temperature has a rather long history [l|^-|33[ . The heavy mesons 
play very important role in our understanding of nonperturbative dynamics of QCD. First determinations of leptonic 
decay constant of pseudoscalar, B c meson at zero temperature were made twenty years ago 0, HE]. Such charged 
meson decays play important role to extract the magnitudes and phases of the Cabbibo-Kobayashi-Maskawa(CKM) 
matrix elements, which can help us understand the origins of CP violation in and beyond the standard model. Our 
aim in this work is to investigate the temperature dependence of mass and leptonic decay constants of the pseudoscalar 
B Cl r] c and rjb mesons taking into account the additional operators arising at finite temperature. The pseudoscalar 
decay constant, fp is defined by vacuum to meson matrix element of the axial vector current as: 



(0|(Q 1 7^75<92) (0)\P)=if P q ) „ (1) 

where Qx,2 = c or b and P — B Cl rj c or r/t- In thermal field theories, the meson mass, mp and its decay constant, fp 
should be replaced by their temperature dependent versions. 

The paper is organized as follows. In section 2, we obtain thermal QCD sum rules for the masses and decay 
constants of the considered pseudoscalar mesons calculating the spectral densities and nonperturbative corrections. 
In section 3, we present our numerical calculations and discussions. 



II. THERMAL QCD SUM RULES FOR DECAY CONSTANTS AND MASSES OF HEAVY 

PSEUDOSCALAR, B c ,r) c AND rjb MESONS 

Taking into account the new aspects of the finite temperature QCD, sum rules for the masses and decay constants 
of the heavy pseudoscalar mesons containing b and/or c quark are derived in this section. The starting point is to 
consider the following responsible two-point thermal correlation function: 

I%T) = i J d 4 xe^ x (T(J P (x)J p \0)) ), (2) 

where T denotes the temperature, T is the time ordering product and J P (x) = Qi{x)i^i^Q 2 {x) is the interpolating 
current of the heavy pseudoscalar mesons. The thermal average of the operator, A = T (J p (x) J p t(0)) appearing in 
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the above correlation function is expressed as: 



(A) = 



Tr(e~' m A) 
Tr{e-P H ) 



(3) 



where H is the QCD Hamiltonian, ft — 1/T is inverse of the temperature T and traces are performed over any 
complete set of states. 

As we previously mentioned, to obtain sum rules for physical observables, we need to calculate the aforementioned 
correlation function in two different ways. In QCD or theoretical side, the correlation function is calculated in deep 
Euclidean region, q 2 <C — ^qcd y i & OPE where the short or perturbative and long distance or nonperturbative 
contributions are separated, 



n QCD {q,T) = Il pert (q 1 T) + Tl nonpert {q 1 T). 



(4) 



The perturbative contribution is calculated using perturbation theory, whereas the nonperturbative contributions are 
expressed in terms of the thermal expectation values of the quark and gluon condensates as well as thermal average 
of the energy density. The perturbative part can be written in terms of a dispersion integral, hence 



U Q CD (q,T) = J 



ds P (s,T) +n nonpert^ T , 

s — q 



(5) 



where, p(s,T) is called the spectral density at finite temperature. The thermal spectral density at fixed |q| is written 



p{q,T) = - ImW ert (q,T) tanh ( ^° 

7T \ Z 



(6) 



In order to calculate the p(q, T) in the lowest order in perturbation theory, we use quark propagator at finite 
temperature [34] as: 



S(q) = (7„ <f + m) 



1 



q 2 — m 2 + is 



+ 2nin(\q \) S(q 2 - to 2 ) 



(7) 



where n(x) — [exp(/3x) + 1] is the Fermi distribution function. Using the above propagator, after some calculations 

(8) 



we find the imaginary part of the correlation function as: 

ImU(q,T) = L(q Q )+L(-q Q 
where, 



L(q ) - -N, 



dk io 2 — k 2 + k • q — LUiqo — TO1TO2 

8^ 



| ^1 - ni(wi)Vl - ti 2 (w 2 )^ + ni(uji)n 2 {u2) S(q - U) X - u 2 ) 
(l - ni(ui))ri2(L)2) + (l - «2(w2))ni(wi) 5(q - ui +w 2 )|. (9) 



Here, u>i 



i\ and lo 2 = y (q-k) 



As it is seen, the L(qo) involves two pieces. The first term, which 



includes delta function 5 (go — wi — ^2) survives at zero temperature and is called the annihilation term. The second 
term, which includes delta function S(qo — W1+W2) is called scattering term and vanishes at T — 0. The delta function, 
5(qo — cui — LU2) in Eq. ^ gives the first branch cut, q 2 > (mi + to 2 ) 2 , which coincides with zero temperature cut 
that describes the standard threshold for particle decays. On the other hand, the delta function, S(qo — ui\ + cj 2 ) 
in Eq. ([9]) shows that an additional branch cut arise at finite temperature, q 2 < (mi — to 2 ) 2 , which corresponds to 
particle absorption from the medium. Taking into account these contributions, the annihilation and scattering parts 
of spectral density in the case, q = can be written as: 



p^(s,T) = P0 (s) [l - n(f (l + ^t)) - n(f ( 



for (mi + m 2 ) 2 < s < 00, and 



s ( m 2 — to 2 
2 V 1 a 



p ^( S) T) = p o(s )[n(f (l + ^t))-n(-f (l 



(10) 



(11) 
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for < s < (mi — TO2) 2 with mi > m2. Here po{s), is the spectral density in the lowest order of perturbation theory 
at zero temperature and is given by: 

M«) = ^q 2 (s)v(s), (12) 

where q(s) = s — (mi — ra-i) 2 and v(s) = \fT— 4mim2/l(s)- 

In our calculations, we also take into account the perturbative two-loop order a s correction to the spectral density. 
For equal quark masses case this correction at zero temperature can be written as (2l| : 

A, (.) = ^ [- - — ( T - 3) + In _ + <?(«)] , (13) 



where -F('u) and G(v) functions have the following forms 

3 

4^ ~ I6v 16 ' 16 



3 21 1Sv 3 ^ 3 



and 



9 1v 2 

G <")-^ + 8-^' < 15 ' 



Here w = v(s) and we replace the strong coupling a s in Eq. (|13[) with its temperature dependent lattice improved 
expression [3, [38|. When doing the numerical calculations for B c meson, the contribution coming from two-loop 
diagrams is used for unequal quark masses case p as [2ll . 23], but since its expression is very lengthy, we do not present 
its explicit expression here. 

To calculate the nonperturbative part in QCD side, we use the nonperturbative part of the quark propagator in an 
external gluon field, A^(x) in the Fock-Schwinger gauge, x^A"(x) — 0. Taking into account one and two gluon lines 
attached to the quark line, the massive quark propagator in momentum space can be written as [2l[ : 



S aa (k) 



— m 4 \k A — m A y 

ft + m 



o- K \(ft + ni) + {ft + m)a K x 



where, 



-g 2 (t c t d ) aa 'G c afi (0)G d ^(0) '"J ^ (/ aftw + f a ^ v + f a ^)(fi + m), 



(16) 



(17) 



In order to proceed, we also need to know the expectation value, {TrG a pG^ u ). The Lorentz covariance at finite 
temperature allows us to write the general structure of this expectation value in the following way: 



1 



(T^GapG^) - —(g a »g/3v - g av gt3n)(G a Xa G aAa ) 
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gafigpv - gavgpy, - 2(u a Upgp l> ~ u a u u gp^ - upu^g al , + upu v g ail ) (u x Q 9 XrT u a ), 



(18) 



where, u M is the four-velocity of the heat bath and it is introduced to restore Lorentz invariance formally in the 
thermal field theory. In the rest frame of the heat bath, u M = (1, 0, 0, 0) and u 2 — 1. Also Q 9 Xcr is the traceless, gluonic 
part of the stress-tensor of the QCD. Therefore, up to terms necessary for our calculations, the non perturbative part 
of massive quark propagator at finite temperature takes the form: 



7 G° 
S {k) - 4 9{t) (k 2 -m 2 ) 2 



cr K \{ft + m) + (ft + rn)a K \ 



i g 2 6 



2 zaa' 



9 (k 2 - m 2 ) 4 

(k 2 - A(k ■ uf \ + (m 2 - 4(fc • uf^j ft + 4(k- u){k 2 ~ m 2 ) i (u a Q 9 a ^)\, 



1 3m(fc 2 + m ft) ( G c^ G ca^ 



(19) 



Using the above expression and after straightforward calculations, the nonperturbative part in QCD side is obtained 
as: 



j-^nonpert 



I! dx i 



(a s G 2 ) 



48tt 



m\{— 1 + x) — [m\ + q 2 (—l + x)^j: 



mf (-1 + x) 4 (5 - 20x + 3x 2 ) + mfm 2 (-l + x) 2 



x (5 - 6x + 12x 2 - 14x 3 + 6x 4 ) - x 4 (^24 m\ q 4 (-l + x) 3 + 6 q 6 (-1 + x) 4 + m\ (-12 + 14a; + 3a; 2 ) 
+ 2 m 4 q 2 (15 - 31a; + 16a; 2 )) + m 4 (-l + x) 2 x( - 2 g 2 (-l + x) 2 (-l + 16a;) + m§(-8 + 58a; - 56a; 2 + 3a; 3 )) 
+ mf (-1 + x)x 2 (24 q 4 (-1 + xfx + m 4 (3 - 45x + 47a; 2 + 3x 3 ) + 3 m\ q 2 (-1 + 23a; - 44.x 2 + 22a; 3 )) 
- m 3 ™ 2 (-l + x)x(q 2 (-8 + 5a; + 9x 2 - 7x 3 + x 4 ) + 2 to 2 (4 - 3x + 9x 2 - 12a; 3 + 6a; 4 ) 
x (3 q 4 (l - 5a; 2 + 6a; 3 - 2x 4 ) - m\ q 2 (6 - 9x 2 + 2x 3 + x 4 ) + m\ (3 + 6a; 2 - 10a; 3 + 6a; 4 )) 



TTll 77J 2 x 



+ 



m 2 (-l +x)- [m\ + q 2 {-l + a;)) : 



72tt 



6 (-1 + ?i 2 )(-l + x) 4 (5 - 20a; + 3a; 2 ) - 6 mf m 2 (-l + u 2 ) 



x (-1 + x) 2 (5 - 6a; + 12x 2 - 14.x 3 + 6a; 4 ) + 2 m\ (-1 + x)x 2 -3 m 4 (-1 + w 2 )(3 - 45a; + 47a; 2 + 3x 3 ) 

- m\q 2 (-1 + x) ( - 9(1 - 22a; + 22a; 2 ) + m 2 (23 - 227a; + 227a; 2 )) + q 4 (-l + xf + x)x + u 2 (-U + 96x 

- 92a; 2 + a; 3 )) - m\{-l + x) 2 a;^6 m^(-l + m 2 )(-8 + 58a; - 56a; 2 + 3a; 3 ) + q 2 (-l + x)(12 (1 - 17a; + 16x 2 ) 
+ m 2 (-37 + 259x - 239a; 2 + 8x 3 ))) + a; 3 6 m\ (-1 + u 2 ) x (-12 + 14a; + 3a; 2 ) + q 6 (-l + x) 3 ( - 36(-l + x)x 
+ u 2 (9- 47a; + 47a; 2 )) + 2 m\ q 4 (-1 + xf ( - 72(-l + x)x + u 2 (9 - 85x + 89x 2 + a; 3 )) + m\ q 2 (-1 + x) 

x (l2(15 - 16a;)x + u 2 (9 - 195x + 215x 2 + 8.x 3 )) +2m\m 2 (-1 + x)x 6 m\ (-1 + m 2 )(4 - 3x + 9x 2 - 12x 3 



+ 6x 4 ) + q 2 (-l + x) ( - 3(8 + 3x - 6x 2 + x 3 ) + 2 u 2 (12 + x - 8x 3 + 4x 4 ) 



2 mi m 2 x A 



3 m\(-l + u 2 ) 



x (3 + 6x 2 - 10x 3 + 6x 4 ) + <7 2 (-l + x) 2 ^9(-l - 2x + 2x 2 ) + m 2 (9 + 14x - 12x 2 - 4x 3 + 2x 4 )) + m\ q 2 (-1 + x) 
x (2(-6 - 6x + 3x 2 + x 3 ) + 2 u 2 {9 + 7x - 8x 3 + 4x 4 )) - 4(-l + x)x|2 m\ m 2 (-1 + x) 2 x (7 - llx + 8x 2 ) 
- m 4 (-l + x) 2 (-25 + 55x - 47x 2 + 8x 3 ) - 2 mi m 2 (-1 + x)x 2 (m 2 (4 - 5x + 8x 2 ) + 2q 2 {-2 + 3x - 2x 2 +x 3 )) 
+ 2 mf (-1 + x) x (ro 2 .(-14 + 29x- 29x 2 ) +q 2 (U- 38x + 44x 2 - 21x 3 + x 4 )) + x 2 ^ 4 (-l + x) 2 (9 - llx 
+ llx 2 ) + m\{9 - 15x + 23x 2 + 8x 3 ) + 2 m\ q 2 (-9 + 22x - 30x 2 + 16x 3 + x 4 )) j (q ■ u) 2 J |, (20) 

where, (G 2 ) = {G c aj3 G cal3 ) . 

Now, we turn our attention to the physical or phenomenological side of the correlation function. The hadronic 
spectral density is expressed by the ground state pseudoscalar meson pole plus the contribution of the higher states 
and continuum. According to quark-hadron duality, the continuum is expected to be well approximated by the QCD 
spectral density calculated in perturbation theory starting at some threshold so- Therefore, the hadronic spectral 
density can be written as: 



o had {s) = Jp 



f 2 P (T)m 4 P (T) 
(mi + mi) 2 



S( S -m 2 P ) + d(s- SQ )pP ert {s) 



(21) 



Matching the phenomenological and QCD sides of the correlation function, sum rules for the mass and decay 
constant of pseudoscalar meson are obtained. To suppress the contribution of the higher states and continuum, the 
Borel transformation over the q 2 as well as continuum subtraction are performed. As a result of the above procedure 
and after lengthy calculations, we obtain the following sum rule for the decay constant: 



fp{T) m 4 P {T) eSF = {m x + m 2 ) 2 



«o(T) 



(m 1 +m 2 ) 2 



ds p a ^ ert {s) e 



(m l -m 2 y 



ds p s - pert {s) e—fr + BIl nonpert y 

(22) 



where M 2 is the Borel mass parameter. 



The sum rule for the mass is obtained applying derivative with respect to — to the both sides of the sum rule 



6 



for the decay constant of the pseudoscalar meson in Eq. (|22|) and dividing by itself: 



m 2 P {T) 



where. 



_> _ f(Sm 2 y ds P a ' pert ( g ) s eM-Th) + fo mi ~' m2rds P s ' pert ( s ) s exp(-^) + n 1 " onpert (M 2 ,r) 



•(mi -m 2 ) 2 



(mi — m 2 ) s 



ds p s -P ert {s) exp(—As) + BTl non P ert 



(23) 



ni nonper *(M 2 ,T) = M 4 -r—r BU nonpert , 
dM z 

and BH nonpert shows the nonperturbative part of QCD side in Borel transformed scheme and is given by: 



(24) 



Bn nonpert = / fa 



96 7T M 6 X 4 (-1 



■ exp 



mfa; — m 2 (— 1 + x) 



{(a s G 2 ) - m\{-l + xf + m\m 2 {-\ + x) A x 



M 2 x(-l + x) 

x (-1 + 2.t)+:e 4 (-12to^ M 4 (-1 + a;) 3 + 12 M 6 (-1 + a;) 4 + 2 m 4 . M 2 x (-1 + a;) - m^a; 2 ) 

+ m\x{~\ + x) 3 (2M 2 {-1 + x) 2 + m 2 (l - 3a; + a; 2 )) + m 2 a; 2 (-l + a;) (l2M 4 a;(-l + x f + m\x(-\ + x + x 2 ) 

+ 3m 2 M 2 (l ~3x + Ax 2 - 2a; 3 )) + mfm 2 a;(-l + x) 2 ( - m\x{\ - 2x) 2 + M 2 (2 - 9x + 6x 2 + a; 3 )) 

- m 1 m 2 (-l + x)x 2 (mjx 2 (l - 2x) - m 2 2 M 2 x{& - 9x + x 2 ) + QM 4 (l+x- 4x 2 + 2a: 3 )) 

+ 3 a s (e 9 ) mf(-l + x) 6 - mfm 2 x(-l + x) A {-l + 2x) + m x m 2 x 3 (-\ + x) (m A 2 x{l - 2x) 

+ 4M 4 (-1 + x) 2 (2 - x + x 2 ) + mj M 2 (-4 + 3x + 5.x 2 - 4a; 3 )) - m 4 x{~l + x) 3 (m 2 (l - 3a; + x 2 ) 

+ 2M 2 (1 -2x + a; 3 )) - m 2 a; 2 (-l + x)(m\x{-l + x + x 2 ) + m|M 2 (5 - 17a; + 24a; 2 - 12a; 3 ) 

+ M 4 (-l + a;) 2 (-l + 15a; - 7a; 2 + 2a: 3 )) + x 3 [mix 3 + M 6 (-l + a;) 3 (9 - 11a; + 11a; 2 ) + 2m 2 M 2 x 

x (-1 + 4a; - 4a; 2 + a; 3 ) - m§M 4 (-l + a;) 2 (-9 + 7a; + a; 2 + 2a; 3 )) + mjm 2 a; 2 (-l + x) 2 

x (m\{\ - 2a;) 2 + M 2 (l + 6a; - 11a; 2 + 4a; 3 )) }, 



(25) 



where, 9 = 0q O . Following [18|, we also use the gluonic part of energy density both obtained from lattice QCD 
[35l [3r| and chiral perturbation theory [37| . In the rest frame of the heat bath, the results of some observables 
calculated using lattice QCD in [3j| are fitted well by the following parametrization for the thermal average of total 
energy density, (0) : 



(9) = 2(0 9 ) = 6 x l(T 6 ea;p[80(T - 0.1)](GeV 4 ), 



(26) 



where temperature T is measured in units of GeV and this parametrization is valid only in the interval 0.1 GeV < 
T < 0.17 GeV. Here, we would like to stress that the total energy density has been calculated for T > in chiral 
perturbation theory, while this quantity has only been obtained for T > 100 MeV in lattice QCD (for details see 
[35l |36|). In low temperature chiral perturbation limit, the thermal average of the energy density is expressed as [37| : 



(6) = (0£) + 3 p, 

where, (0£J) is trace of the total energy momentum tensor and p is pressure. These quantities are given by: 



\ 1*1 270 Fi \T 



p = 3 T 



771^ T 
2 7T 



15 T 105 T 2 



+ 



8 m„ 128 ml 



exp 



m^ 
T 



(27) 



(28) 



where A p = 0.275 GeV, F n = 0.093 GeV and m n = 0.14 GeV. 

In our calculation we use the temperature dependent continuum threshold, sq(T), gluon condensate, (G 2 ) and 
strong coupling constant, as presented in [lj| (for details see [l3l. liH l35l. [36l l38j). 
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III. NUMERICAL ANALYSIS 

In this section, we numerically analysis the sum rules for the masses and decay constants of the heavy-heavy 
pseudoscalar mesons. We use the values, m c = (1.3 ± 0.05)GeV, nib = (4.7 ± 0.1)GeV and (0 | \ol s G 2 | 0) = 
(0.012 ± 0.00A)GeV' 1 for quark masses and gluon condensate at zero temperature. From the sum rules for the masses 
and decay constants it is clear that they also contain two auxiliary parameters, namely continuum threshold, so and 
Borel mass parameter, M 2 as the main inputs. These are not physical quantities, hence the physical observables should 
be independent of these parameters. Therefore, we should look for working regions for these parameters at which 
the dependence of the masses and decay constants on these parameters is weak. The continuum threshold, sq is not 
completely arbitrary, but it is in correlation with the energy of the first exited state with the same quantum numbers 
as the considered interpolating currents. We choose the values 44 GeV 2 < s < 46 GeV 2 , 11 GeV 2 < s < 12 GeV 2 
and 94 GeV 2 < so < 97 GeV 2 for the continuum threshold in accordance with B c , rj c and rjb channels, respectively. 
The working region for the Borel mass parameter, M 2 is determined as following. Its lower limit is calculated requiring 
that the higher states and continuum contributions constitute approximately 30% of the total dispersion integral. Its 
upper limit is obtained demanding that the mass sum rules should be convergent, i.e., contribution of the operators 
with higher dimensions is small. As a result of the above procedure, the working region for the Borel parameter is 
found to be 10 GeV 2 < M 2 < 25 GeV 2 , 6 GeV 2 < M 2 < 12 GeV 2 and 15 GeV 2 < M 2 < 30 GeV 2 in B c , rj c and rj h 
channels, respectively. 

Our calculations show that in the working regions the dependence of the considered observables on auxiliary pa- 
rameters is weak. We depict the dependence of masses and decay constants on the temperature, T in Figs. [THf)] These 
figures contain the results obtained using both lattice QCD and chiral perturbation parametrization for the gluonic 
part of the energy density. These figures depict that both parametrization of lattice QCD and chiral perturbation 
theory predict the same result in validation limit of lattice QCD fit parametrization, i.e., 0.10 GeV <T< 0.17 GeV. 
These figures also show that the masses and decay constants remain unchanged approximately up to T ~ 100 MeV, 
but after this point, they start to diminish with increasing the temperature. Near the critical or deconfincmcnt tem- 
perature, the decay constants reach approximately to 38% of their values in vacuum, while the masses are decreased 
about 5%, 10%, 2% comparing with their values at zero temperature for B c , rj c , rjb mesons, respectively. From these 
figures, we obtain the results on the decay constants and masses at zero temperature as presented in Tables I and II. 
The quoted errors in these Tables are due to the errors in variation of the continuum threshold at zero temperature, 
Borel mass parameter as well as errors coming from fit parametrization of the temperature dependent continuum 
threshold, gluon condensate and strong coupling constant and uncertainties existing in other input parameters. These 
Tables also include the existing predictions of the other works as well as experimental data. The Table II depicts a 
very good consistency between our results and the experimental data on masses but from Table I, we see that the 
present work results and the results existing in the literature (see Table I) on the decay constant are comparable up 
to presented errors. 

Our results for the leptonic decay constants at zero temperature as well as the behavior of the masses and decay 
constants of the considered pseudoscalar heavy mesons with respect to the temperature can be checked in the future 
experiments. The obtained behavior of the observables in terms of temperature can be used in analysis of the results 
of the heavy ion collision experiments. 





f Bc (MeV) 


fvc(MeV) 




Present Work 


476 ± 27 


421 ± 35 


586 ± 61 


QCD sum rules [19, 24, 251 


400 ± 25 


350 




Potential Model [28] 


400 ± 45 


402 


599 


Lattice QCD Method [30] 


489 ±7 


438 ± 11 


801 ± 12 


Experiment [39] 




335 ± 75 





TABLE I. Values of the leptonic decay constants of the heavy-heavy pseudoscalar, B c , ij c and rjb mesons in vacuum. These 
results have been obtained using the values M 2 = 15 GeV 2 , M 2 = 6 GeV 2 and M 2 = 20 GeV 2 ior B c , rj c and rj b particles, 
respectively. 
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m Bc (GeV) 


m„ c (GeV) 


m Vb (GeV) 


Present Work 


6.37 ±0.05 


2.99 ±0.04 


9.58 ± 0.03 


Experiment [40J 


6.277 ± 0.006 


2.9803 ±0.0012 


9.3909 ± 0.0028 



TABLE II. Values of the mass of the heavy-heavy pseudoscalar, B c , r\ c and mesons in vacuum. The same values as Table I 
for the auxiliary parameters have been used. 
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FIG. 1. The dependence of the mass of B c meson on temperature for Chiral and Lattice QCD parametrization of the gluonic 
part of the energy density. 
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FIG. 2. The same as Fig. [Qbut for f Br 
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FIG. 3. The same as Fig. [T]but for m Vc . 
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